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Abstract 


Iron is commonly used as a structural and shielding material in nuclear devices. 
The accuracy of its nuclear data is critical for the design of nuclear devices. The 
evaluation data of “Fe isotopes in the latest version of the CENDL-3.2 library from 
China was significantly updated. This new data must be tested before it can be used. To 
test the reliability of this data and assess the shielding effect, a shielding benchmark 
experiment was conducted with natural Fe spherical samples using a pulsed deuterium- 
trittum neutron source at the China Institute of Atomic Energy (CIAE). The leakage 
neutron spectra from the natural spherical iron samples with different thicknesses (4.5, 
7.5, and 12 cm) were measured between 0.8-16 MeV after interacting with 14 MeV 
neutrons using the time-of-flight method. The simulation results were obtained by 
Monte Carlo simulations by employing the Fe data from the CENDL-3.2, ENDF/B- 
VIILO, and JEDNL-5.0 libraries. The measured and simulated leakage neutron spectra 
and penetration rates were compared, demonstrating that the CENDL-3.2 library 
performs sufficiently overall. The simulation results of the other two libraries were 
underestimated for scattering at the continuum energy level. 

Keywords: Iron sphere, CENDL-3.2, shielding benchmark experiment, pulsed 14 
MeV neutron source, time-of-flight method. 


1. Introduction 


A complete set of nuclear data is frequently required as an input for simulation 
calculations for the design of reactors, nuclear fuel cycles, nuclear medicine, and 
applications of nuclear technology!!! Most international countries with nuclear power 
have created evaluation libraries to meet the application requirements of nuclear data. 
The United States, Western Europe, Soviet Union, Japan, International Atomic Energy 
Agency, and China have constructed nuclear data libraries since the 1960s, such as the 
ENDF/B, JEEF, BROND, JENDL, FENDL, and CENDL, respectively. The latest 
version of the CENDL-3.2 library was released in 2020"!, whereas that of the JENDL- 


5.0 library was released in 20211, 

Nuclear evaluation data refers to the results obtained from assessments or 
evaluations of nuclear processes or events and can be used to assess the safety and 
efficacy of nuclear technology. This type of data may include the measurements of 
radioactive materials, radiation levels, and other parameters related to nuclear events or 
processes. These data are widely used for microscopic cross-sectional measurements!” 
71, as well as for accelerator and reactor designs!*!, Nuclear evaluation data are important 
for assessing nuclear safety, particularly for developing and operating nuclear power 
plants and other facilities that handle radioactive materials. Nuclear technologies 
require careful monitoring and assessment to ensure that safety risks are minimized and 
the technologies are used effectively and responsibly!!. A benchmark experiment is 
critical for determining the accuracy of the assessment data and is an important aspect 
of analyzing nuclear data before it is distributed to users!'®!, Macroscopic sample 
measurements are an important prerequisite for benchmark experiments; for certain 
critical evaluation data, a series of macroscopic sample measurements must be 
performed to verify their reliability. Benchmark experiments are important for 
validating the nuclear data, which is essential for a wide range of nuclear applications. 
A benchmark experiment involves comparing predictions made using nuclear data with 
experimental measurements obtained under controlled conditions. By comparing the 
results of these calculations with the experimental data, scientists can evaluate the 
accuracy of the nuclear data in predicting the behavior of nuclei under various 
conditions. 

This process is critical for ensuring the safety and reliability of nuclear 
applications because inaccurate nuclear data can lead to errors in the reactor design and 
operation, as well as present potential safety risks!!! Therefore, benchmark 
experiments are critical for improving the understanding of nuclear processes and 
developing more accurate models for predicting the behavior of nuclei. 

Iron plays a critical role in the nuclear industry owing to its unique properties!'7!, 
including its ability to absorb and decelerate neutrons, which are fundamental particles 
used in nuclear reactions. Iron is extensively used in the construction of nuclear reactors, 
particularly in the form of steel, an alloy mainly composed of iron. In nuclear reactors, 
iron-based materials are used to manufacture various components, such as reactor 
pressure vessels, steam generators, and fuel cladding. These components are crucial for 
the safe operation of reactors and are critical for controlling nuclear reactions. Iron is 
also used for nuclear fuel production. In certain nuclear reactors, iron is used as a matrix 
material to encapsulate the nuclear fuel pellets. Encapsulation helps contain radioactive 
material and prevents it from leaking into the environment. Overall, the role of iron in 
the nuclear industry is essential because it provides a reliable and durable material for 
constructing nuclear facilities and manufacturing nuclear components. Improving the 
accuracy of the nuclear evaluation data of Fe isotopes, especially °°Fe, is highly 
significant for reducing the redundancy of the shielding design of nuclear devices and 
improving their miniaturization and refinement. The penetration problem, which is the 
movement of particles in a wide geometric region packed with matter, is one of the 
most common challenges in radiation shielding calculations. To improve the accuracy 


of the calculation results of radiation shielding after the neutrons deeply penetrate the 
iron, high-quality nuclear evaluation data of Fe is required as supporting information. 

Several shielding benchmark experiments of Fe have been conducted to test the 
evaluation data of the Fe isotopes. These include measurements of the leakage neutron 
spectra of the Fe sphere by a D-T neutron source in the Lawrence Livermore National 
Laboratory (LLNL)!"1, the angular leakage neutron spectra of slab samples in the Japan 
Atomic Energy Research Institute (JAERI)""*!, and the leakage neutron and y spectra of 
slab samples in the Technical University of Dresden (TUD)! 1. The absolute reaction 
rates in spherical Fe systems were measured at the Institute of Nuclear Physics and 
Chemistry of CAEP''®'7! and the leakage neutron spectra of the slab samples were 
measured at the China Institute of Atomic Energy (CIAE)""®!, These results provide 
important experimental data for analyzing and improving the Fe evaluation data in their 
libraries”, 

To assess the quality of the “Fe evaluation data in the CENDL-3.2 libraries and to 
provide high-quality evaluation data for the deep penetration shielding calculations of 
Fe, a shielding benchmark experiment of Fe spheres based on a fusion D-T neutron 
source was conducted in this study. Using the time-of-flight method, the leakage 
neutron spectra of D-T neutrons passing through iron spheres of different thicknesses 
were measured, which were 4.5, 7.5, and 12-cm-thick, and the measured energy ranged 
between 0.8-16 MeV. The MCNP-4C program was used to obtain the simulated results; 
the evaluation data of Fe from the CENDL-3.2, ENDF/B-VIII.0, and JEDNL-5.0 
libraries were used. The simulated and experimental results were compared and the 
differences were analyzed. The results demonstrate that the simulated results of the 
CENDL-3.2 library are in sufficient agreement with the experimental results, and the 
evaluated nuclear data of Fe from the CENDL-3.2 library can be used as the optimal 
data for nuclear device shielding calculations. 


2. Experiment Setup 


The Chinese Academy of Atomic Energy Sciences established an integral 
experimental system for verifying the nuclear evaluation data], Because there is 
only one reaction channel for neutron and hydrogen (H) elastic scattering, and its cross- 
section is commonly used as the international standard, the leakage neutron spectrum 
obtained by measuring the hydrogen-containing material (standard sample) is the main 
source for checking the integration of the experimental system. The CIAE system was 
experimentally tested by measuring the leakage neutron spectra of polyethylene and 
water samples'*>!, The schematic of the experiment for measuring the leakage neutron 
spectrum of an Fe sphere sample is shown in Figure 1, which mainly includes a neutron 
source, sample, collimated shield, detectors, and an electronic acquisition system. The 
primary experimental parameters are listed in Table 1. 
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Fig.1 Overall experimental layout. 
Table 1. Main experimental parameters 
Parameter Description 
Neutron source 14 MeV D-T neutron, Neutron yield 2x10°/s 
Inner Radius 4.85+0.2cm, External radius 
9+0.1cm, 9368+2g 
Inner Radius 4.5+0.1cm, External radius 
9+0.1cm, 20.14+0.2 kg 
Inner Radius 9+0.1cm, External radius 
16.5+0.1cm, 166.12+0.2 kg 


Polyethylene sphere 
Fe sphere No.1 


Fe sphere No.2 


Distance between 
neutron source and 378.4+0.2 cm 
pre-collimator 
Distance between 
neutron source and 770.27+0.2 cm 
detector 


2.1 Neutron source 

Pulsed 300 keV D ions were generated using a Cockcroft-Walton accelerator at 
the Key Laboratory of Nuclear Data in CIAE, which bombards a tritium-loaded 
titanium (T-Ti) target to produce pulsed neutrons of approximately 14 MeV through the 
D(T, n)*He reaction'**!, The diameter of the T-Ti target was 2.2 cm, and that of the active 
zone was 1.6 cm; the thickness was approximately 2.6 um. The parameters related to 
the D-T neutron source are as follows. The beam intensity, pulse width, frequency, and 
neutron yield are approximately 10 uA, 2.5 ns, 1.5 MHz, and 2 x 10°/s, respectively. 
The neutron yields in this experiment were acquired using the associated particle 
method, which measured the “He particles using an Au-Si surface barrier detector 
situated 135° from the beamline 27731, 
2.2 Samples 

Polyethylene and natural iron were selected as the samples; among these, the 
polyethylene sample was used to test the accuracy of the data received from the 
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measurement system. The spherical Fe shell model of the sample is shown in Figure 2, 
which is joined by two hemispherical shells. There are two holes in the spherical shell, 
which are used to accommodate the target tube of the neutron source; the tube is used 
to hold the associated *He particles to enable the T-Ti target to penetrate the center of 
the shell. The dimensions are listed in Table 1, and the mass ratios of the various 
elements of the sample iron sphere are listed in Table 2. A new spherical shell (Fe sphere 
No. 3) is used in these experiments, which is composed of Fe spheres Nos.1 and 2. 
Considering the influence of the two hemispheres on the leakage neutron spectrum at 
the fitting point, the fitting plane is rotated by approximately 15° when the sample is 
placed. Thus, the detector is a complete block of iron spheres at a stereo angle to the 
sample. 


Unit: cm 


Fe sphere No. 1 Fe sphere No.2 


Fig. 2 Model of spherical iron shell. 
Table 2. Composition of the sample 


Nuclide Mass ratio 
natFe 99.9% 

C 0.015% 

Si 0.03% 
Mn 0.02% 

Al 0.02% 
Pb+Sn+Bi+Ca+S+P 0.015% 


2.3 Collimator 

During the experiment, neutrons not only interact with the sample, but also with 
other devices in the laboratory; these neutrons may enter the detector after scattering 
and form background counts. To reduce the influence of the background, a system that 
combines collimating and shielding was used in the experiment, which was composed 
of a pre-collimator and an additional collimator embedded in the wall°1. By using this 
system and placing the detector outside of the collimator in the wall, the background 
spectrum measured in the experiment was far lower than the other spectra, as shown in 
Figure 3. 
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Fig. 3 A comparison of the results after the normalization of each sample. 


2.4 Detector systems and data processing 

The time-of-flight spectra of the leaking neutrons were measured after passing 
through the spherical Fe sample. An EJ-301 liquid scintillator detector that was 5.08 
cm in diameter and 2.54 cm thick was used; the distance between the detector and 
neutron source was 76941 cm at the 90° direction of the beam. Figure 4 presents the 


electronic circuit used in the experiment. 
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Fig. 4 Diagram of the electronics circuit. 
The acquisition system using the CAMAC bus unified the EJ-301 detector and 
beam pick-up signals into the KMAX system to obtain the charge integration of the 
signals (QDC), pulse shape discrimination spectrum (PSD), and time-of-flight 


spectrum (TOF) of the leakage particles. 
The QDC of the signal was used as the energy (En) of the particle, after which the 


energy scale was calibrated with '°’Cs and 7*Na y radioactive sources. 

Considering PSD, particle discrimination is performed according to the different 
times of the pulse crossing the zero-point generated by the different particles; the 2160A 
module produced by CANBERRA accomplished this function. 

The TOF is derived from the anode signal of the detector as the starting time signal. 
There is a copper ring in the accelerator tube that generates an induced charge when a 
D-ion pulse beam cluster passes through it to form a pick-up signal, which is used as 
the stop signal after an appropriate delay. Therefore, to determine the actual time of the 
neutron flight, the time of one cycle is deducted from the observed time. 

The neutron and y were screened and selected by an offline analysis using the 
Kmax software to obtain the neutron and y time-of-flight spectra of the EJ-301 detector, 
where the y time-of-flight spectra can be used to obtain the time distribution of the 
pulsed neutrons. 

During the experimental data processing, the final time-of-flight spectrum was 
obtained using the following equation: 

TOF Final — TOFinXEEX = TOF ouT (1) 
KXNaXS NnXS 

Here, TOF;, is the actual neutron spectrum entering the detector, €y is the detector 
efficiency, K is the coefficient when measuring the neutron yield by the accompanying 
particle method, Ng is the number of alpha particles obtained from the measurement, S 
is the detector area, TO Foyr is the neutron event detected by the detector, and N, is the 
count of the neutron source. The spectrum of the final time-of-flight was obtained as 
the neutron flux per unit neutron per unit area (cm?). 


3. Simulation 


In addition to the experimental measurements, experimental simulations are also 
an important part of the shielding benchmark research. The simulation was 
implemented by using the MCNP code. The neutron time-of-flight spectra simulated 
results were obtained using the neutron evaluation data of a specific library and MCNP 
code, which can be directly compared to the experimental results. Subsequently, the 
simulated results were compared with the experimental results to evaluate and verify 
the evaluation data. The neutron source, sample, and detector are the most important 
parameters in the experiment. An accurate description of the neutron source (energy 
spectrum distribution, angular flux distribution, and pulsed time distribution) and 
detector (neutron detection efficiency curve) is critical for ensuring the accuracy of the 
simulation results), 

3.1 Energy spectra and angular distribution of the neutron source 

MCNP simulations require input information of the neutron source; the precision 
of the description of the neutron source directly affects the accuracy of the output results. 
In general, the energy spectrum and angular distribution of the neutron source can be 
obtained using both experimental measurements and simulations. However, accurate 
experimental measurements are difficult given the situation of an accelerator in an 
experimental hall; therefore, the TARGET program was used to calculate the energy 
spectrum of the neutron source?'!, The TARGET program, developed by PTB in 


Germany, can calculate the energy spectrum and angular distribution of a deuterium— 
tritium neutron source based on the incident D-particle and target parameters (energy, 
pulse width, and tritium target thickness). The results of the TARGET calculations are 
shown in Figure 5, where the neutron energy decreased as the exit angle increased. The 
spectrum of the neutron source and the angular distribution calculated by using the 
TARGET program can be directly used as the parameters of the neutron source in the 


MCNP program. 
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Fig. 5 The energy spectrum and angular distribution of the D-T neutron source. 


3.2 Time distribution of the pulsed neutron 


In the simulation of the time-of-flight spectrum, the pulsed time distribution of the 
neutron source requires a precise description. The pulse-time distributions of the 
neutron sources for the different samples were obtained by extracting the time-of-flight 
spectra of the leakage y-particle from the measurements. In the offline analysis of the 
data, the y-particle was selected to extract the time-of-flight spectrum, which can be 
used as the time distribution of the pulsed neutrons after appropriate processing, as 


shown in Figure 6. 
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Fig. 6 Time distributions of pulsed neutrons of the different samples. 

3.3 Detector efficiency 

A liquid scintillator detector was used to measure the energy spectrum of the fast 
neutron in the experiments. The detector efficiency curve is an important parameter for 
measuring the neutron spectrum by liquid scintillator detectors. The detector efficiency 
curve was calculated using the NEFF program developed by PTB (Germany)!*”!, which 
was obtained after describing the detector geometry, light response curve, and energy 
threshold in the NEFF program. The threshold energy for simulating the detector 
efficiency curve was chosen to be 0.3 times that of the Compton edge energy of !°’Cs, 
and the neutron energy corresponding to this Compton edge was 0.8 MeV. The 
simulation results are shown in Figure 7, which were stored as the detector efficiency 
in the input card of the MCNP code. 
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Fig.7 Simulation result of the detector efficiency. 

3.4 The time-of-flight spectrum 

The energy spectra, angular flux distribution, pulse-time distribution of the source 
neutrons, and detector efficiency curve of the liquid scintillator detector are described 
in the input card of the MCNP-4C codel*l. The environment of the experimental hall 
was complicated and difficult to describe for the simulated physical model; thus, only 
neutrons within the solid angle of the detector to the sample were considered. According 
to the background measurement shown in Figure 3, the neutrons scattering into the 
detector by the walls, ground, and other facilities can be ignored. Structures that 
significantly influence the neutrons of the incoming detector include the target, sample, 
and collimators. Therefore, the entire experimental hall was simplified as an air-filled 
cylinder. The target, collimators, and sample were modeled according to the actual 
structure and dimensions. The data for polyethylene and Fe were adopted from the 
evaluation data of the CENDL-3.2, ENDF/B-VIIL0"*!, and JEDNL-5.0 libraries, and 
all other structural material information of the nuclei were obtained from the evaluation 
data of the CENDL-3.2 library. The simulated results of the time-of-flight spectra can 
be directly obtained using the MCNP program. A ring detector was used in the 
simulation, and the time-of-flight spectrum obtained from the simulation indicates the 
neutron flux per unit area (cm°), which can be directly compared with the 
experimentally measured spectrum. The calculation model is shown in Figure 8. 
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Fig.8 Simulation model in the MCNP program. 


4. Uncertainties Analysis 


The final time-of-flight spectrum obtained in the experiment is given by Equation 
1, where the coefficient k of the associated particle method is corrected by comparing 
the measurements and simulations of the neutron source; thus, the value of K is as 
follows: 

NEx 
J NsimX€smXNanX5S (2) 
where Ngy is the neutron count obtained by measuring the time-of-flight spectrum of 
the neutron source, Ngn is the corresponding associated a count, €s,y is efficiency of the 
detector obtained from the simulation, Ns;yis the flux of the simulated time-of-flight 
spectrum of the neutron source, and S is the detector area. According to Equations 1 
and 2, TOFinai is as follows: 
TOFfinat = e st (3) 

Because the experimental data are finally normalized using the relative coefficients, 
considering that the ratio of the experimental measurements of the neutron source to 
the simulation results is used as the relative coefficient K, most of the systematic errors 
are reduced or eliminated, including the angular errors caused by the measurements of 
the associated a particle and time-of-flight spectra. Only the statistical error from the 
simulation of the neutron source (<0.5%), the statistical uncertainties in the 
measurement of the time-of-flight spectrum of the neutron source (<0.5%), the 
statistical uncertainties in the measurement of the alpha particles (<0.2%), the relative 
detection efficiency uncertainties (<3%), and the statistical uncertainties in the neutron 
counts of each channel remain in the final uncertainties >), 


5. Experimental results and discussion 


5.1 Neutron source measurement results 
First, the TOF spectrum of the D-T neutron source was measured directly in the 
experiment, the results of which are shown in Figure 9(a). The leakage neutron 


spectrum is composed of the following two parts: the direct contribution of the neutrons 
emitted by the source and the contribution of the scattered neutrons. The direct 
contribution from the neutron source was mostly between 120 and 206 ns, and the 
contribution from the scattered neutrons was mostly between 206 and 620 ns. The target 
tube-scattered, air-scattered, and collimator-scattered neutrons accounted for the 
majority of the scattered neutrons. 
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Fig.9 (a) Comparison of the simulated and experimental results of the time-of-flight spectra of the 

neutron source. (b) Comparison of the simulated and experimental results for the polyethylene samples. 

The model shown in Figure 8 was used to simulate the time-of-flight spectrum of 
the neutron source. Table 3 lists the integral values of the experimental and simulated 
results for the two intervals and the C/E ratio (simulated value/experimental value). The 
experimental and simulated results for the direct contribution interval of the neutron 
source were significantly similar, indicating that the measurement system was generally 
reliable. The experimental results are higher than the simulated results for the scattered 
neutron interval; Figure 9(a) demonstrates that there is a peak near 350 ns 
(approximately 2.5 MeV), which is primarily owing to the fact that the unreacted 
deuterium (D) ions will be continuously deposited in the target during the experiment, 
and subsequently react with the incident D ions behind with the D(d, n)*He reaction to 
produce neutrons of approximately 2.5 MeV. The probability of this reaction increases 
as the time of the target usage increases. In addition, the situation in the experimental 
hall is more complex than that in the simulation model. It is impossible to consider all 
scattered neutrons, leading to lower simulation results. 

Table 3. Comparison of the experimental measurements and simulations of the neutron sources 


Proportion in Proportion in 


Interval Simulation ; à Experiment i C/E 
simulation experiment 
Neutron source 1.07E-08 94.70% 1.07E-08 92.30% 1.000+0.030 


Scattered neutrons 6.01E-10 5.30% 8.93E-10 7.70% 0.673+0.020 


5.2 Measurement results of the polyethylene sample 

To further ensure the reliability of the measurement system, measurements of the 
leakage neutron spectrum were obtained, and the simulation of a standard sample 
(polyethylene sphere) was performed. Data from the CENDL-3.2, ENDF/B-VIII.0, and 
JEDNL-5.0 libraries were used for the polyethylene material in the simulation; 


simulation and experimental results of the different libraries were obtained, as shown 
in Figure 9(b). 

By calculating the energy distribution of the emitted neutrons after their interaction 
with C and H nuclides, the spectrum of the leaking neutrons after the interaction of the 
14 MeV neutrons with the polyethylene spheres can be approximately divided into three 
intervals: 

(1) The elastic scattering interval (120-200 ns, depending on the pulse width of 
the source neutron) primarily includes C and H elastic scattering. The minimum energy 
of a neutron emitted by the elastic scattering of C is approximately 10.5 MeV (180°). 

(2) The vast majority of the H elastic scattering interval (200-250 ns) is owing to 
H elastic scattering. 

(3) The low-energy neutron interval includes inelastic scattering from C, the C 
reaction (n, n'3a), elastic scattering from H, and multiple scatterings. 

Table 4 lists the integral values of the experimental results for each interval, the 
integral values of the simulation results for different libraries, and the C/E values. 

Table 4. Measurement results of the polyethylene sample 


Interval Integral Value RES ee C/E 
rate 
Elastic scattering interval 6.84E-09 60.93% 
; H elastic scattering interval 8.77E-10 7.81% 
Experiment | 
low-energy neutron interval 1.59E-09 14.16% 
0.8-16MeV 9.31E-09 82.93% 
Elastic scattering interval 6.48E-09 57.74% 0.947+0.029 
i ing i 8.62E-10 7.68% 0.983+0.031 
CENDL-32 H elastic scattering HER o 
low-energy neutron interval 1.48E-09 13.18% 0.93 140.029 
0.8-16MeV 8.82E-09 78.60% 0.948+0.029 
Elastic scattering interval 6.57E-09 58.57% 0.961+0.030 
ENDF/B- H elastic scattering interval 8.77E-10 7.81% 1.000+0.031 
VIILO low-energy neutron interval 1.44E-09 12.80% 0.905+0.028 
0.8-16MeV 8.89E-09 79.18% 0.955+0.029 
Elastic scattering interval 6.48E-09 57.73% 0.947+0.030 
i ing i 8.64E-10 7.70% 0.985+0.031 
JEDNL-5.0 H elastic scattering ameeval o 
low-energy neutron interval 1.51E-09 13.46% 0.951+0.030 
0.8-16MeV 8.86E-09 78.89% 0.951+0.029 


The cross-section of the hydrogen elastic scattering was used as the standard. A 
comparison of the results shown in the preceding table indicates that the deviation of 
the C/E values of CENDL-3.2 library and JENDL-5.0 in the H elastic scattering interval 
is within 3%, and the deviation of the simulation result from ENDF/B-VIII.0 is less 
than 1%, confirming the reliability of the measurement system. Notably, the simulation 
results obtained by using the ENDF/B-VIIL.0 library in the low-energy neutron interval 
are significantly underestimated by approximately 9% when compared to the 
experimental results, which is primarily because the neutron energy spectrum of the (n, 
n'3a) reaction of C obtained by the ENDF/B-VIII.0 library is significantly lower than 


that of the other libraries. 
5.3 Measurement results of the Fe sample 

The time-of-flight spectra of the leaking neutrons from the deuterium-tritium 
pulsed neutrons interacting with the spherical Fe samples of three different thicknesses 
(4.5, 7.5, and 12 cm) were measured using the time-of-flight method. The MCNP 
program was used to calculate the simulated results of the leaking neutron spectra of 
the spherical Fe samples using the following three libraries: CENDL-3.2, ENDF/B-VIII, 
and JEDNL-5.0. Figure 10 presents a comparison of the experimental and simulated 
results. 
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Fig.10 Comparison of the simulated and experimental results of iron spheres with varying thicknesses: 
(a) 4.5 cm, (b) 7.5 cm, and (c) 12 cm. 

The 14 MeV neutrons interacted with the natural iron sample through four main 
reaction channels: elastic scattering (n, el), discrete energy level inelastic scattering (n, 
inl)D, continuous energy level inelastic scattering (n, inl)C, and (n, 2n) reactions. The 
NDPlot software can be used to draw the secondary neutron energy spectrum of a 
natural iron sample incident with the 14 MeV neutrons, as shown in Figure 11. By 
combining the pulse time widths that were experimentally obtained, the time-of-flight 
spectrum of the outgoing neutron was divided into four intervals based on the energy 
intervals of the different reaction channels, as listed in Table 5. 


—natFe Total 
— natFe (n,el) 
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Fig.11 Secondary neutron energy spectrum of the natural iron samples at the 14 MeV neutron incidence 
given by the CENDL-3.2 library. 

According to Figure 11, the entire experimental energy spectrum is divided into 
the following four intervals: 12.5-16 MeV indicates elastic scattering (n,el), 8.5—12.5 
MeV indicates discrete energy level inelastic scattering (n,inl)D, 2.4-8.5 MeV indicates 
continuous energy level inelastic scattering (n,inl)C, and 0.8—2.4 MeV indicates the 
(n,2n) reaction. 


Table 5. The energy interval corresponding to each reaction channel 


Reaction Energy Interval Time Interval 
(n,el) 12.5 MeV~16 MeV 139 ns~157 ns 
(n,inl)D 8.5 MeV~12.5 MeV 157 ns~190 ns 
(n,inl)C 2.4 MeV~8.5 MeV 190 ns~359 ns 
(n,2n) 0.8 MeV~2.4 MeV 359 ns~622 ns 


The experimental and simulated spectra were integrated according to the 
aforementioned energy intervals; the experimental integrated values, simulated 
integrated values, C/E values, and neutron penetration rates with different thicknesses 
were obtained for each interval, as listed in Table 6. 

Table 6. C/E values and penetration rates of the iron spheres with different thicknesses 


Thickness Interval Experiment Penetration rate 


(nel) 5.25E-09 46.78% 
(n,inl)D 1.16E-09 10.35% 
45cm  (ninl)C 2.06E-09 18.36% 
(n,2n) 3.08E-09 27.42% 
Total(0.8-16MeV) _1.16E-08 102.92% 
(nel) 3.62E-09 32.25% 
(n,inl)D 1.08E-09 9.60% 
75cm (n,inl)C 1.05E-09 9.35% 
(n,2n) 1.87E-09 16.68% 
Total(0.8-16MeV)  7.62E-09 67.87% 
(nel) 2.02E-09 17.98% 
(n,inl)D 6.50E-10 5.79% 
12cm (n,inl)C 1.06E-09 9.47% 
(n,2n) 2.33E-09 20.77% 
Total(0.8-16MeV 6.06E-09 54.00% 


Table 7. C/E values of the different thicknesses 


Simulation results C/E 
Fhickness Interval ENDF/B- ENDF/B- 

CENDL-3.2 VILO JENDL-5.0 CENDL-3.2 VILO JENDL-5.0 

(n,el) 5.34E-09 5.32E-09 5.34E-09 1.018+0.031 1.01340.031 1.017+0.031 

(n,inl)D 1.14E-09 1.14E-09 1.15E-09 0.979+0.031 0.984+0.031  0.988+0.031 

45cm  (n,inl)C 2.10E-09 1.82E-09 1.92E-09 1.021+0.032  0.882+0.028 0.934+0.029 
(n,2n) 3.19E-09 3.13E-09 3.13E-09 1.036+0.032  1.016+0.032  1.017+0.032 
Total(0.8-16MeV) 1.18E-08 1.14E-08 1.15E-08 1.019+0.031 0.988+0.031  0.999+0.031 

(n,el) 3.52E-09 3.52E-09 3.55E-09 0.973+0.030 0.972+0.03 0.980+0.030 

(n,inl)D 1.09E-09 1.08E-09 1.10E-09 1.015+0.032  1.007+0.032  1.022+0.032 

7.5cm  (n,inl)C 1.06E-09 9.07E-10 9.63E-10 1.015+0.032 0.864+0.028  0.918+0.029 
(n,2n) 2.00E-09 1.89E-09 1.92E-09 1.067+0.033  1.008+0.031  1.025+0.032 
Total(0.8-16MeV) 7.68E-09 7.40E-09 7.53E-09 1.008+0.031  0.971+0.03 0.988+0.031 

(n,el) 1.88E-09 1.86E-09 1.89E-09 0.930+0.029  0.921+0.029 0.934+0.029 

(n,inl)D 7.92E-10 7.88E-10 8.04E-10 1.220+0.038 1.214+0.039  1.237+0.039 

12cm  (n,inl)C 1.09E-09 8.98E-10 9.66E-10 1.029+0.032 0.845+0.027  0.909+0.029 
(n,2n) 2.53E-09 2.34E-09 2.40E-09 1.086+0.034  1.002+0.031  1.031+0.032 
Total(0.8-16MeV) 6.30E-09 5.88E-09 6.06E-09 1.039+0.032  0.970+0.030 0.999+0.031 


The following can be concluded based on the results shown in Tables 6 and 7: 

(1) The neutron penetration rate in the 0.8-16 MeV interval slightly exceeds 1 at a 
thickness of 4.5 cm, as shown in Table 6. This is because the thickness is relatively 
small, the neutrons have not decelerated to below 0.8 MeV, and the neutrons present a 
certain proliferation owing to the (n,2n) reaction. As the thickness of the sample 
increases, the neutrons constantly decelerate, and the neutrons below 0.8 MeV become 
more abundant, while the neutrons in the 0.8-16 MeV interval continue to decrease, as 
well as the penetration. 

(2) The neutron scattering interval ranges between 4-16 MeV, which includes 


elastic scattering, discrete energy-level inelastic scattering, and continuous energy-level 
inelastic scattering. Because the leaking neutrons in this part of the interval are 
primarily affected by the cross-section of the scattering neutrons, the transmittance 
exhibits an exponentially decaying trend. The exponential decay law satisfied by the 
penetration rate should be y = e74 , where y is the penetration rate and x is the 
thickness of the Fe sample. The fitting results using the experimental data in Table 6 
are shown in Figure 12. According to the results, the fitting parameter A is 0.15495 for 
the elastic scattering energy region and 0.08423 for the energy region between 2.4—-16 
MeV. 
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Fig.12 Trend of the neutron penetration rate in different energy regions with varying thicknesses of the 
spherical iron. 

(3) The variation in the C/E values of the 0.8-16 MeV neutrons with different 
spherical iron thicknesses using the various libraries is shown in Figure 13. The total 
cross-section provided by the CENDL-3.2 library is apparently slightly larger, and that 
provided by the ENDF/B-VIII.0 library is smaller; this deviation increases as the 
thickness increases. The simulation results of the JENDL-5.0 database are in sufficient 
agreement with the experimental results, and the overall deviation can be controlled 
within 3%. 
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Fig.13 The C/E value trends with varying iron thicknesses in the interval between 0.8-16 MeV. 

(4) Figure 3 and Table 6 demonstrate that when the thickness increases from 4.5 
cm to 12.5 cm, the leakage neutrons in the elastic scattering interval gradually decrease. 
When the thickness is increased from 4.5 cm to 7.5 cm, the leakage neutron count in 
the inelastic scattering interval significantly reduces; however, when the thickness is 
increased from 7.5 cm to 12 cm, it does not significantly change. This is owing to the 
fact that the amount of leakage neutrons in the inelastic scattering interval increases as 
the high-energy neutrons decelerate and decreases when the thickness increases; the 
change in the neutron count as the thickness changes is a result of both effects. These 
two effects also apply to the counting law for leaking neutrons in the (n,2n) interval. 
Moreover, only neutrons with energies greater than 0.8 MeV are measured in this 
experiment, thus, a large fraction of neutrons are decelerated to below 0.8 MeV and are 
not included in the count. The overall trend of the leakage neutron counts in the (n,2n) 
interval decreased and then increased as the thickness increased. 

(5) The C/E values obtained by these three libraries gradually decreased as the 
thickness increased, whereas the C/E values for the discrete energy-level inelastic 
scattering intervals gradually increased. This suggests that the cross-sections of the 
inelastic scattering for the data from all three libraries may be slightly larger, and the 
cross-sections of the elastic scattering are slightly smaller. 

(6) The simulated values of the inelastic scattering interval at the continuum 
energy level for both the JENDL-5.0 and ENDF/B-VIIL.0 libraries were approximately 
10% lower than the experimental values for all three thicknesses. This indicates that the 
cross-sections of the inelastic scattering at the continuum energy level are small, which 
are given by the JENDL-5.0 and ENDF/B-VIIL.0 libraries for the 14.5 MeV incident 
neutron natural iron samples. As shown in Figure 14, the energy spectrum of the 
secondary neutrons given by the ENDF/B-VIII.0 library is larger between 0-3 MeV but 


smaller between 3-8 MeV, which may be the reason for the underestimation of the 
simulation results using the ENDF/B-VIILO library in this interval. The simulation 
results of the CENDL-3.2 library exhibited different degrees of bias in the simulated 
values in the (n,2n) interval. 
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Fig.14 Secondary neutron spectrum of the inelastic scattering at the continuum energy level produced 


by 14 MeV incident neutrons on natural Fe. 


6. Summary 


In this study, an integration experiment of a sphere Fe shield was conducted using 
a D-T pulsed neutron source at the China Institute of Atomic Energy. The leakage 
neutron spectra of Fe spheres with thicknesses of 4.5, 7.5, and 12 cm were measured 
using the time-of-flight method. To test the reliability of the iron data in the CENDL- 
3.2, ENDF/ B-VILO, and JENDL-5.0 libraries, the time-of-flight spectra were simulated 
using the MCNP-4C code. In the simulation, in addition to the physical model, the 
detector efficiency, energy and angular distributions of the neutron source, and pulse 
neutron time distributions are fully considered, such that the time-of-flight spectrum 
demonstrates an improved fine structure, which will help compare the time-of-flight 
spectra for short time intervals. The entire time-of-flight spectrum was divided into four 
intervals: elastic scattering (n, el), discrete energy level inelastic scattering (n, inl)D, 
continuous energy level inelastic scattering (n, inl)C, and (n, 2n) reaction intervals. 
Subsequently, the simulation and experimental results at different intervals were 
compared to verify the data for the different reaction channels. 

According to the comparison results, the following conclusions can be drawn: (1) 
The simulation results of the CENDL-3.2 library are in the best agreement with the 
experimental results. (2) The simulations from the JEDNL-5.0 library were 
approximately 8% lower in the continuum energy level inelastic scattering interval, 


whereas the simulations from the ENDF/B-VIII.0 library were approximately 10% 
lower. This may be owing to the fact that the ENDF/B-VHI.0 and JENDL-5.0 libraries 
underestimate the energy spectrum of the secondary neutrons of the high energy part of 
the continuum energy level inelastic scattering. 
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